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ABSTRACT 

A rare coincidence of scales in standard particle physics is needed to explain 
why A or the negative pressure of cosmological dark energy (DE) coincides with 
the positive pressure Pq of random motion of dark matter (DM) in bright galaxies. 
Recently Zlosnik et al. (2007) propose to modify the Einsteinian curvature by 
adding a non-linear pressure from a medium flowing with a four- velocity vector 
field If^. We propose to check whether a smooth extension of GR with a simple 
kinetic Lagrangian of If^ can be constructed, and whether the pressure can bend 
space-time sufficiently to replace the roles of DE, Cold DM and heavy neutrinos in 
explaining anomalous accelerations at all scales. As a specific proof of concept we 
find a Vector-for-A model (VA-model) and its variants. With essentially no free 
parameters, these appear broadly consistent with the solar system, gravitational 
potentials in dwarf spiral galaxies and the bullet cluster of galaxies, early universe 
with infiation, structure formation and BBN, and late acceleration with a 1:3 ratio 
of DM:DE. 

Subject headings: Dark Matter; Cosmology; Gravitation 

The incompleteness of standard physics and Einstein's General Relativity (GR) is ev- 
ident from the smallness of the cosmological constant A or the vacuum energy density 
~ (O.OOleV)^, compared to the expected quantum pressure ~ (lO^^eV)"^ at scales 

of the Planck mass nip = \J^- Current speculations of the new physics of A are as free 
as analogous speculations of the Pioneer Anomaly (Turyshev, Nieto, Anderson 2006), both 
represent acceleration discrepancies of order ~ 7ao, driven by unidentified (likely unrelated) 
pressures ~ 72Po, where Oq = 1.2 A/ sec^, and Pq = are scales of acceleration and pres- 
sure. On intermediate scales, galaxy clusters and spiral galaxies often reveal a discrepant 
acceleration of order (0.1 — 2)ao. GR, if sourced primarily by baryons and photons with 
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negligible mass density of neutrinos and other particles in the Standard Model or variations, 
appears an adequate and beautiful theory in the inner solar system, but appears increasingly 
inadequate in accounting for astronomical observations as we move up in scales from lOOAU 
to 1 kpc to IGpc. The universe made of known material of positive pressure should show 
a de-accelerating expansion as an open universe, but instead it is turning into accelerating 
now, evidenced by much dimmer supervonae detected at redshift unity. A standard rem- 
edy to restore harmony with GR and fit successfully large scale observations (Spergel et al. 
2006 and references therein) is to introduce a "dark sector", in which two exotic compo- 
nents dominate the matter-energy budget of the Universe at the redshift z with a split of 
VLbe '■ ^DM = 3 : (1 + z)^ approximately: a Dark Energy (DE) as a negative pressure and 
nearly homogeneous field described by unknown physics, and a Cold Dark Matter (DM) as a 
colissionless and pressureless fluid motivated by perhaps the MSSM physics. However, antic- 
ipating several new particles from the LHC, the success of this Concordance Model still gives 
little clue to the physics of governing the present 1 : 3 ratio of its constituents. This ratio 
is widely considered improbable, because standard particle physics expects a ratio 1 : 10^^°. 
Here we speculate whether the 3 : (1 + ratio could come from a coincidence of scales of 
Oq = I.2A/ sec^ with the cosmological baryon energy density p^c^ ~ 3.5 x (1 + z)^Pq. 

A deeper link of DM and DE: It is curious that the distribution of DM in dwarf galaxies 
is extremely ordered, something that the cuspy ACDM halos are still struggling to explain 
even with the maximum baryonic feedback (Gnedin & Zhao 2002). For example, on galaxy 
scales the Newtonian gravity of DM = ~ Qb and Newtonian gravity of baryons 
93 = have a tight correlation: 

{qVobT - g"" ^a"^, 9 = 9dm + 93. (1) 

where n > 1 (Zhao & Famaey 2006). This rule holds approximately at all radii R of all spiral 
galaxies of baryonic mass Mb{R) and circular velocity Vc{R) within the uncertainty of the 
stellar mass-to-light and object distance. For low surface brightness galaxies or at the very 
outer edge of bright spirals, the gravity g is weaker than ao, our empirical formula predicts 
9"^ /9b = (y/^ / RY / {GMb/ R^) = V^/{GMb) ~ ao, which is essentially the normalisation 
of (baryonic) TuUy-Fisher relation (McGaugh 2005). Bulges and central part of elliptical 
galaxies are dominated by baryons inside a transition radius where the baryon and DM 
contribute about equally to the rotation curve, Eq.([T]) predicts goM = 9b = ao/2; we 
can define a DM pressure Pq = x at transition by multiplying the local gravity 
{9dm + 9b) = ^0 with the DM column density above this radius = ao/(87rG). This scale 
Pq appears on larger scales too. All X-ray clusters have gas pressure and the DM random 
energy density comparable to Pq. The amplitude of the scale ao appears in the r~^ cusp of 
CDM halos too (Xu, Wu, Zhao 2007, Kaplinhat & Turner 2001). These can be understood 
since the last scattering shell at z = 1000 has a thickness 2L ~ lOMpc and contains typical 
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potential wells of depth (? jN ~ (1000 km s"^)^ due to inflation, where N = 10^, hence 
the typical internal acceleration is c^/N/L ~ 0.2ao. Also a DM sphere of radius 5Mpc 
turning non-linear now would fall in with an acceleration ~ 200 x Hq x 5Mpc ~ Qq. While 
correlations of baryon and DM can generally be understood in a galaxy formation theory 
where DM and baryons interact, the unlimited freedoms of dark particles means a good 
spread of its concentration, hence the correlation would have substantial history-dependent 
variance from galaxies to galaxies and radii to radii. For example, DM is unexpected in Tidal 
Dwarf Galaxies, is observed for its Oq acceleration (Gentile et al. 2007). The tightness of 
such hidden regulations on DM at all radii for all galaxies is anomalous, at least challenging 
in the standard framework. 

It is even more curious that DM in various systems and DE are tuned to a common scale 
Pq, hence requiring a coincidence in two dark sectors. These empirical facts are unlikely 
random coincidences of the fundamental parameters of the dark sectors. Since all these 
anomalies are based on the gravitational acceleration of ordinary matter in GR, one wonders 
if the dark sectors are not just a sign of an overlooked possible field in the gravitational 
sector. 

Continue along Zhao (2006), here we propose to investigate whether the roles of both 
DM and DE could be replaced by a vector field in a modified metric theory. This follows 
from two long lines of investigations pursued by Kostelecky, Jacobson, Lim and others on 
consequences of symmetry-breaking in string theory, and by Milgrom, Bekenstein, Sanders, 
Skordis and others driven by astronomical needs. These two independent lines were first 
merged by the pioneering work of Zlosnik et al.(2007). Existence of an explicit Lagrangian 
satisfying main constraints for the solar system, galaxy rotation curves and cosmological 
concordance ratio remains to be demonstrated. 

Warming up to vector field: In Einstein's theory of gravity, the slightly bent metrics for 
a galaxy in an uniform expanding background set by the fiat FRW cosmology is given by 

g^^dx^dx" = -(1 + -^)d{cty + (1 - —)a{tfdf (2) 

where dt^ — [dx^ + dy^ ■\- dz'^) is the Euclidian distance in cartesian coordinates. In the 
collapsed region of galaxies, the metric is quasi-static with the potential ^{t, x, y, z) — 
^{t, X, y, z) due to DM plus baryon, which all follow the geodesies of g^^,. 

Modified gravity theories are often inspired to preserve the Weak Equivalence Principle, 
i.e., particles or small objects still go on geodesies of above physical metric independent of 
their chemical composition. Unlike in Einstein's theory, the Strong Equivalence Principle 
and CPT can be violated by, e.g., creating a preferred frame using a vector field. The 
Einstein-Aether theory of Jacobson & Mattingly (2001) is such a simple construction, where 
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a unit vector field is designed to couple only to the metric but not matter directly. 
It has a kinetic Lagrangian with linear superposition of quadratic co-variant derivatives 
V(c^f/)V(c^f/), where c^U^ is constrained to be a time-like four-momentum vector per unit 
mass by —g^JJ^U" = 1. The norm condition means the vector field introduces up to 3 new 
degrees of freedom; e.g., a perturbation in the FRW metric (Eqj2]) has c^U^ = g^yC^U'' ~ 
(c^ -|- containing a four- vector made of an electric-like potential $ and three 

new magnetic-like potentials. But for spin-0 mode perturbations with a wavenumber vector 
k, we can approximate — (1, 0) ~ which contains just one degree of freedom, i.e., 

the flow potential V{t,x,y,z). We expect an initial fluctuation of c|k|y ~ |$| ~ c^N^^ = 
10~^c^ can be sourced by a standard infiaton; the vector field tracks the spectrum of metric 
perturbation (Lim 2004). 

Most recently Zlosnik et al. (2007) suggested to replace the linear AVf/Vf/ with a non- 
linear kinetic Lagrangian F(AVf/Vt/) to extend Jacobson's framework. They showed this 
class of non-linear models is promising to produce the DE effect in cosmology and the DM- 
like effect in the weak field limit. Here we continue along the lines of the pioneering authors, 
but aim for a single Lagrangian with parameters in good match with basic observations of a 
range of scales. 

A Simple Lagrangian for A: The difficulty of writing down a specific Lagrangian is that 
there are infinite ways to form pressure-like terms quadratic to co- variant derivatives of the 
vector field. Simplicity is the guide when choosing gravity since GR plus ACDM largely 
works. Let's start with forming two pressure terms for any four- momentum-like field A^^ 
with a positive norm mc^ = —gapA'^AI^ by 

8.GJ(.4) . 1 ' , S.GKiA) . Yl^Yl^ (3) 

3 \ m J mm 

where the RHSs are co-variant with dimension of acceleration squared, and Vy = A'^Va or 
Vq, stands for the co-variant derivative with space-time coordinates along the direction of 
the vector A or the dummy index a respectively. From these we can generate two simpler 
pressure terms K and J of the unit vector field U" by 

J = J{U) ~ 0, K = KiU) ~ ^ in galaxies 
~ %i!g ' ~ universe 

where the approximations hold for with negligible spatial components and nearly fiat 
metric (EqE]). Note the J and K are constructed so that we can control time-like Hubble 
expansion and space-like galaxy dynamics separately^ The i^-term, with a characteristic 
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pressure scale = Pq in galaxies, is the key for our model. The J-term, meaning critical 
density, has a characteristic scale N'^Pq ~ IO^'^Pq: at the epoch of recombination z — 1000 
when baryons, neutrinos, and photons contribute ~ (8, 3, 5) x IO^Pq respectively to the term 
J = ^fjg ; so the epochs of equality and recombination nearly coincide. 

Now we are ready to construct our total action S — J d'^x] —g\ in physical coordinates, 
where the Langrangian density 

C = + Lr,, + Lj + LK+ {UM''+ l)L^, (5) 

where R is the Ricci scalar, is the ordinary matter Lagrangian. For the vector field part, 
L"^ is the Lagrangian multiplier for the unit norm and we propose the new Lagrangian 




where the non-negative continuous functions Xi{x) = [0, X{x) — X{N)]^^^^ , X{x) = (l + |) 
where the subscript i = either n or A^. Incidentally, n = gives GR. The cutoffs (e.g., with 
n = ±1) guarantee a bounded Hamiltonian with kinetic terms Lk and Lj always bounded 
between ±N^Po (e.g., in a lab near Earth K - (lO^^ - 10'^^)Po > N^Po, so Lk = 0). The 
condition at the tidal boundary K — J — is well-behaved too (cf. eq. 44-48 of Famaey 
et al. 2007 on Cauchy problem). Note 1 - ^ > i^rnin = (1 + iV/n)"" ~ 10"^^ and 
1 - > A^s = (1 + N/N)-'' ~ 2-3. 

Taking variations of the action with respect to the metric and the vector field, we can 
derive the modified Einstein's equation (EE) and the dynamical equation for the vector 
field. The expressions are generally tedious (Halle 2007), but the results simplifies in the 
perturbation and matter-dominated regime that interest us. As anticipated in Lim (2004) 
the ij-cross-term of EE yields ^ — $ = for all our models, which means incidentally twice 
as much deflection for light rays as in Newtonian. As anticipated in Dodelson & Liguri 
(2007), the ti-term EE can be casted into that of an unstable harmonic oscillator equation 
with a negative string constant V + biHV — (1 — iiB)b2H^V = S{^, ^) if (1 — /^b) > 0, so 
we expect that HV tracks $. The tt-term of the EE takes the form 

SnGp = SplbH^ + 2V • [(1 - A„) V$] - Aq - g($, V, V) (7) 



SwGKia = 2g"''\c^yaU~<){c^Vi3U^) - 2{c^V aU'^){c^V 13U") which change details of structure formation, 
PPN, and gravitational waves, which are beyond our goal here. 
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where we approximated 1 — Xn{x) ~ 2~" = /i^ as a constant in matter dominated regime 
where J < N'^Pq and the Q-term is zero for static galaxies and uniform FRW flat cosmology. 
So the tt-equation of Einstein reduces to the simple form 



AnGp = V^<l>-V 



A„ ( 1 V'l' 

ao 



in galaxies 



= if 2 _ matter-dominated FRW (9) 

Here the pressure from the vector field creates new sources for the curvature. The term 
^^'^"JnG^'^^ Poisson equation acts as if adding DM for quasi-static galaxies. A cosmo- 

logical constant in the Hubble equation is created by 

^"^^ fhM) . (10) 



SwG ' > " ' („ _ l)(„ - 2) 

For binary stars and the solar system, AirGp — ^ is true because the gravity at 
distances 0.3AU to 30AU from a Sun-like star is much greater than the maximum vector 
field gradient strength Nao, so ^ = 0; in fact, |V$| ^ ^ ~ (10^ - 10^)ao, and the 
typical anomalous acceleration is NaoPmin ~ 10~^°ao, well-below the current detection limit 
of lO^'^ao (Soreno & Jezter 2006). This might explain why most tests of non-GR effects 
around binary pulsars, black holes and in the solar system yield negative results; Pluto at 40 
AU and the Pioneer satellites at 100 AU might show interesting effects. Extrapolating the 
analysis of Foster & Jacobson (2006), we expect GR-like PPN parameters and gravitational 
wave speeds in the inner solar system. 

Near the edges of galaxies, we recover the non-relativistic theory of Bekenstein & Mil- 
grom (1984) with a function 

p{x) = 1 - Xn{x) ~ Pmin + X, if X = ^ < 1. (11) 

Note that fi{x) — x hence rotation curves are asymptotically flat except for a negligible 
correction fimin ~ 10^^^. In the intermediate regime a; = 1 our function with 1 — A„(a;) ~ 
(0.55 — 0.6) for n = 2 — 5 respectively. Eq.([T]) argues that galaxy rotation curves prefer 
a relatively sharper transition than /i(x) = x/{l + x) = 0.5 at x = 1 (Famaey, Gentile, 
Bruneton, Zhao 2007) where we can identify gB/igoM + Qb) = p{x). So our model should 
fit observed rotation curves. 

A 2 

For the Hubble expansion: the vector field creates cosmological constant-like term 
9Pq below the zero-point of the energy density in the solar system because the zero point 
of our Lagrangian (Eql6]) is chosen at N'^Pq < K < +oo. During matter domination, the 
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contribution of matter SnGp and Aq to the Hubble expansion (Eq|9]) is further scaled- up 
because the effective Gravitational Constant Ge// = G/fis = 2"^ > G with GR being the 
n = special case. B Coming back to the original issue of the 3 : 1 ratio of matter density to 
our cosmological constant, Eq.dni) predicts that : ~ — : Mllfli^ which is close 

to the desired 3 : (1 + 2;)^ ratio. Adding neutrinos makes the explanation slightly poorer. 
So the DE scale is traced back to a separate coincidence of scale, i.e., the present baryon 
energy density p^c^ ~ APq, where Pq contains a scale for the anomalous accelerations on 
galactic scale. Our model predicts that DE is due to a constant of vacuum, preset by the 
modification parameter n of the gravity; n = gives GR. 

In our model, the effective DM (the dog) follows the baryons (the tail) throughout the 
universal (1 + z)^ expansion with a ratio set by n. To fit the ACDM-like expansion exactly, 
we note the Hubble equation for a flat FRW cosmology with vector field and standard mix of 
baryons, neutrinos and photons ~ ^'m^^ ~ o 000025 ~ -'■ yields at the present epoch 



1 - = nt^^'^ (12) 



The 2nd equality fixes /i^^ = 2" = (8 — 8.4) if we adopt ao/c ^ Hq/6 ^ 12km/s/Mpc and 
fl^^'^ = (0.25 — 0.3). The 1st equality would predict an uncertain but very small neutrino 
mass m-i, ~ ±0.3eV. 

The BBN also anchors any modification to GR. In the radiation-dominated era \J\ = 
^fj^ 3> N'^Po, the dynamics is driven by 

SnGp ~ — Aq — Atv, in radiation-dominated FRW, (13) 

where = — j^N{x)d{PQN'^x^) = —N'^Pq/8 for n = 3 is a finite negative number, much 
smaller than the radiation pressure ~ (z/1000)^A^^Po- So the early universe is GR-like, 
especially the Hubble parameter at the BBN, insensitive to the precise value of N'^Po- 

Note a more general version of our Vector-for-A model has a Lagrangian 

Lk + Lj = XMUxi) + XjJiUXj) - PoV{Xk, Aj), (14) 

1 1 

with four vector degrees in UXJ and the scalar Xk/Xj; it is optional to replace UX^ with 

1 

UX^ to reduce the total freedom to 4 as in Bekenstein's TeVeS. Our simple model is equiva- 
lent to the special case of two non-dynamical scalar fields Xk and Xj with ;^ ~ ;^ — ^ 0, hence 



^Whilc Carroll & Lim (2004) found a scale-down of G because they were interested in stable spin-0 modes 
with (1 — /xb) < for a restricted class (04 = 7^ ci) of Jacobson's models, Dodelson & Liguri (2006) argue 
that an unstable growth of the vector field is helpful to structure growth in many gravity theories. 
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K, = }C{U) = K and J = J{U) = J (Eql3]). The potential is smooth with PoV{\k, Aj) = 
iLn [^(Z^-- + A;, - A) - ^n{X -Xj- 2-")] PndX, where P„ = (A"^ - ifn^Po and n{y) 
is the Heaviside function of y. A vector field ^ (mc^ +m$, mcA) with a mass scale m has 
a quantum degeneracy pressure limit ~ ^^"^^ It is intriguing that our model suggests the 
existence of a zero point vacuum energy ~ -PoV(l, 1) ~ 9Po ~ (O.OOleV)^. And the (pos- 
itive) radiation pressure at the epoch of baryon-radiation equality coincides with the cutoff 
energy density PoV(0, 0) ~ —N'^Pq ~ — (0.3eV)^, and the vacuum-to-cutoff energy density 
ratio ~ ~ 10~^ coincides with the cosmic baryon-to-photon or baryon-to-neutrino 

number ratio 77 ~ 3 x 10~^° due to a tiny asymmetry with antibaryons. Can theories like 
quantum gravity and inflation explain these coincidences'? Understanding these might give 
clues to how the four-vector potential of photons decouples from the baryon current vector, 
and decouples from our E&M-like vector field A'^ in spontaneous symmetry breaking in the 
string theory (Kostelecky & Samuel 1989, Carroll & Shu 2006, Ferreira et al. 2007). 

Massive neutrinos are optional for our model because the Lj term creates a massive 
neutrino-like effect in cosmology without affecting galaxy rotation curves. There are a few 
ways to create the impression of a fluid of 2eV neutrino in clusters of galaxies as well (Angus 
et al. 2007, Sanders 2005, Zlosnik et al. 2007). E.g., a general Lagrangian with A/" ~ n would 
have new dynamical freedoms /i = 1 — Ax and 1 — Aj, which satisfy second order differential 
equations in time in galaxies, reminiscent of fluid equations for DM. Then the Bekenstein- 
Milgrom /i-function would acquire a history-dependent non-local relativistic correction of 
order j^f^ ~ 1 if the temporal variation (relaxation) time scale r of the scalar field Xk is 
comparable to the Hubble time. This dynamical correction is hard to simulate, but is most 
important at the tidal boundary of (merging) systems where a condensate of the dynamical 
freedoms A^ and Aj oscillate rapidly, could in principle act as an extra DM source to explain 
some outliers to the Bekenstein-Milgrom theory, e.g., the merging bullet cluster with its 
efficient lensing and high speed (Angus & McGaugh 2007). A dynamical field A j is desirable 
as an infiaton to seed perturbations (Kanno & Soda 2004). 

In summary, we demonstrate as a proof of concept that at least one alternative La- 
grangian for the gravity (Eq JSfT^ can be sketched down to resemble the OR plus ACDM 
but with somewhat less-fining in terms of fitting several types of observations from dwarf 
spiral galaxies to the cosmic acceleration. The keys are a zero-point pressure scale Pq at the 
edge of galaxies, and a universal convergence source term ^f^ic'^^aU^f below the cutoff 
pressure N'^Pq, which is near the epoch of equality and the last scattering. However, the 
CMB should be sensitive to the fin = 2~" modification parameter, [fl It should be feasible 
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to falsify the present model and variations by simultaneous fits to the supernovae distances 
and the CMB. 
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